Hall Effect of Light 
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We derive the semiclassical equation of motion for the wave-packet of light taking into account 
the Berry curvature in momentum space. This equation naturally describes the interplay between 
orbital and spin angular- momenta, i.e., the conservation of total angular- momentum of light. This 
leads to the shift of wave-packet motion perpendicular to the gradient of dielectric constant, i.e., 
the polarization-dependent Hall effect of light. An enhancement of this effect in photonic crystals 
is also proposed. 
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The similarity between geometrical optics and parti- 
cle dynamics has been the guiding principle to develop 
the quantum mechanics in its early stage. One can con- 
sider a trajectory or ray of light at the length scale much 
larger than the wavelength A. By setting h = c = 1, the 
equation for the eikonal ip in the geometrical optics reads 
[<9 M V( X )] 2 = [V</>(r)] 2 — n(r) 2 = 0, where — (r,t) is 
the four-dimensional coordinates, ip(x) = 4>(r) — n(r)t, 
and n(r) is the refractive index slowly varying within 
the wavelength This equation is identical to the 

Hamilton- Jacobi equation for the particle motion by re- 
placing ip by the action S. The equation for the ray of 
light can be derived from the eikonal equation in par- 
allel to the Hamiltonian equation of motion. Deviation 
from this geometrical optics is usually treated in terms 
of the diffraction theory. However, most of the analyses 
have been done in terms of the scalar diffraction the- 
ory, neglecting its vector nature. However, the light has 
the degree of freedom of the polarization, which is repre- 
sented by the spin S — 1 parallel or anti-parallel to the 
wavcvcctor k. It is known that this spin produces the 
Berry phase when the light is guided by the optical fiber 
with torsion Q ■ Also the adiabatic change of the polar- 
ization even without the change of k produces the phase 
change called Pancharatnam phase |3|. In this Letter, 
we will show that the trajectory or ray of light itself is 
affected by the Berry phase which leads to various 
nontrivial effects including the Hall effect. 

For simplicity, we focus on an isotropic, nonmagnetic 
medium; the refractive index n{r) is real and scalar, but 
is not spatially uniform. Let us consider the wavefunc- 
tion for the wave-packet with the wavevector centered at 
k c and position centered at r c . Because of the conjugate 
relation of the position and wavevector, both of them 
inevitably have finite width of distribution. Therefore, 
the relative phase of the wavefunctions at k and k + dk 
matters, which is the Berry connection or gauge field A& 
defined by [A&]aa' = -ie\ k Vke\'k, where e\k is the po- 
larization vector of A-polarized photon, and A = ± corre- 
spond to the right- and left-circular polarization. Here, 



it is noted that the SU(2) gauge field A& is 2 x 2 matrix. 
The corresponding Berry curvature (BC) or field strength 
is given by fifc = x A& + iAk x A&. Due to the mass- 
lessness of a photon, it is diagonal in the basis of the right 
and left circular polarization as given by flk = <r^k/k 3 
[a massive case will be discussed later] . It corresponds to 
the field radiated from the monopole with strength ±1 
located at k = 0. We pick up a correction to the geo- 
metric optics up to the linear order in | Vri(r c )|/fc c . An 
equation of motion (EOM) is derived including an effect 
of Afc in the similar way as the Bloch waves of electrons 
5j. By considering the effective Lagrangian for the cen- 
ter coordinates r c and kf^the variational principle leads 
to the following EOM 



k 

r c = v(r c )-^ + k c x (z c \fl kc \z c ), 

k c = -[Vv(r c )]k c , 
z c ) = -ik c ■ Afejz c ), 



(1) 

(2) 
(3) 



where v(r) = l/n(r) is the velocity of light, and \z) = 
t [z + . Z-] represents the polarization state. This is the 
central result of this Letter, from which the Hall effect of 
light is derived below in a straightforward way. 

The limit of geometrical optics is obtained by neglect- 
ing the terms containing A& or $7^, reproducing Fermat's 
principle. Let us study phenomena caused by A& and 
f2fc. First, the A& term in Eq. J2J causes the phase shift 
by the directional change of the propagation discussed in 
ref. Q, as seen in the following way. The equation for 



\z c ) gives the solution |z° ut ) = [e 
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where 



\z?) 



is the initial state of polarization. O 



is the solid angle made by the trajectory of momentum: 
= § dk ■ [K k ]++ = J s dSk ■ [flk]++ where dSk is the 
surface element in fc-space and S is a surface surrounded 
by the trajectory. This is the phase shift by the direc- 
tional change of propagation 0]. 

On the other hand, the second term in Eq. Q induces 
a new effect found in this Letter; the trajectory of light 
is affected by flk- We note that this term guarantees 
the conservation of total angular-momentum of photon 
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jz = [ r c x fc c + (z c |<73|z c )|^] z , assuming the rotational 
symmetry around the z-axis. From Eqs. Q-Q, one can 
prove dj z /dt = when e(r) = x 2 + y 2 , z). Note that 
the orbital angular-momentum r c x k c is defined only 
when the position r c is well-defined, which necessarily 
leads to the finite distribution of the wavevector and 
hence to the Berry phase. This argument applies to a 
light beam, whereas it does not for the plane wave state. 

This Berry-phase term induces the Hall effect of light. 
Let us first consider reflection and refraction of light at a 
flat interface between two media with different dielectric 
constants. Let z = be the boundary between the two 
dielectric media with n = vq/v± where vq and v\ are 
the velocities of light in the media below and above the 
boundary, and r c and k c of the incident light are within 
the y = plane (Fig. P). 




FIG. 1: Transverse shift of light beams in the refraction and 
reflection at an interface. 

The SnelPs law tells us smOr/vi = shvOi/vq for the 
transmitted light and Or = tt — 9i for the reflected light, 
where 9i,t,r are the angle between the z-axis and k c of 
the incident, transmitted and reflected lights. This rela- 
tion can be derived from Fermat's principle in geometri- 
cal optics. The question is "what is the deviation from 
this Snell's law due to the Berry phase ?" The answer is 
the transverse shifts of the centers of the reflected and the 
transmitted wave-packets. This phenomenon could be re- 
garded as the "Hall effect of light" , wherein the "force" 
to the light, i.e., the change of the dielectric constant, 
is perpendicular to the shift. The shift of the transmit- 
ted light can be easily calculated from the EOM. Strictly 
speaking, the EOM can be applied only when the re- 
fractive index n varies much slower than the wavelength; 
thus the EOM is not reliable when the interface is sharp. 
Nonetheless, it turns out to give the correct answer for 
the shift of the transmitted light even for the sharp in- 
terface, because the shift is governed by the conservation 
law of the ^-component j z of total angular-momentum. 
We assume that j z is conserved in reflection and in refrac- 
tion separately: j\ = j'J, j* = j R , where the superscripts 
/, T and R stand for incident, transmitted and reflected 



lights, respectively. This is required when the light is 
regarded as a quantized object, photon; each photon is 
reflected or refracted with some probabilities. From this 
conservation of j z , the transverse shifts of transmitted 
and reflected wave-packets are estimated as 

A [(z^\a 3 \z^) cos 9 A - (z I c \a 3 \z I c )cos9 I ] 

^ (0/) = kJ^Fj ' (4) 

where A = T or R, k l c is the momentum of the inci- 
dent light, \z I c ,T,R ) stand for the polarizations of the in- 
cident, transmitted and reflected photons. We note that 
\z c ) is obtained from \z c ) and the transmittances (re- 
flectances) of two eigenmodes with linear polarizations. 

It has been already known that the totally reflected 
light beam undergoes shifts of position in the plane of 
incidence (longitudinal) or normal to the plane of inci- 
dence (transverse). The longitudinal shift was first stud- 
ied by Goos and Hanchen \f\. The transverse shift was 
predicted by Fedorov J8J , and have been observed exper- 
imentally by Imbert I sl, fo llowed by a number of theo- 
retical approaches [lOL lllj . The transverse shift occurs 
not only in total reflection, but also in partial reflection 
and refraction [l^. While there are several papers on 
theoretical predictions [l2l ITil Hif , they contain rather 
complicated calculations and the issue remains still con- 
troversial. The shift of the transmitted light has not been 
experimentally measured to the authors' knowledge. Our 
derivation based on the EOM is much clearer and sim- 
pler. Also the result is supported by the numerics de- 
scribed below, showing a wide applicability of our EOM. 

In order to test our prediction Eq. (JIJ , we have numer- 
ically solved the Maxwell equation. It is found that both 
the reflected and transmitted wave-packets are shifted 
along the y or — y direction depending on the circular po- 
larization, which perfectly agrees with Eq. Q as shown 
in Fig.|3 This phenomenon is analogous to the spin Hall 
effect jl5llf 6| . where the Berry phase from band structure 
gives rise to spin-dependent motion of electrons. 

It is noted that the similar effect is expected also for 
massive particle waves with intrinsic angular-momentum. 
However, the effect is suppressed by the mass factor, to. 
For example, for the spin-1/2 massive field, the BC in 
helicity basis is given by tt k = ^ [ (a x e$ + cr 2 e c / ) ) + 

<T3e/c], where E k = Vk 2 + m 2 and ei~,o,<j> are the unit 
vectors of the spherical coordinate system in fc-space. 

The BC of light is proportional to 1/fc 2 . Consequently, 
the transverse shift is a fraction of the wavelength, and 
is small for the visible light. In the following we propose 
a method to enhance the Hall effect of light by use of 
photonic crystals 01- ^ s i s known, the periodic mod- 
ulation of the dielectric constant leads to formation of 
Bloch waves of light analogous to those of electron in 
solids. The degeneracy of right- and left-circular ly po- 
larized light is generally lifted. Hence, the SU(2) gauge 
invariance breaks down to U(l) (Abelian). This is similar 
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(a) n = v 1 V, = 2.0 

X\~2 ln x 50 C 
X ,. ~ 2 x 25 /„. 



(b) n = i' / v. 
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FIG. 2: Transverse shift in the center of the transmitted and 
reflected wave-packets of light for (a) n = 2.0 and (b) n = 0.8 
with the incident wave-packet being right-circularly polarized. 
The solid and dashed lines for the analytic results Eq. fi). 
The filled circles and squares are the results of simulation. 
Xl = 2-K/kl is the wavelength of incident wave-packet. l un i t 
is an arbitrary length-scale. The sign of the shifts is reversed 
for the left-circularly polarized light. 



to the anomalous Hall effect in ferromagnetic materials 
|l8|. In photonic crystals, monopoles in fc-space can exist 
away from k = 0. The anomalous velocity is enhanced 
when k is near the monopole. 




(b) Band structure = 0.5) 
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FIG. 3: (a) Dielectric function and (b) band structure of the 
2D photonic crystal. The Brillouin zone is shown in Figs. [I] 
and|3b) 

For an explicit example, we consider a two-dimensional 
(2D) photonic crystal with the dielectric function e(r), 
■ELi[(e-cos(6 l .r + f)) 2 + (e + 



1 _ 4 

t(r) 3(5+12151+8^) 



008(6* • r - f ))], where b, = (^, -g), b 2 = (0, ff ), 
and 63 = — bi — b 2 . The spatial distribution of e(r) is 
shown in Fig.^a). For simplicity, we shall focus on the 
case with k z = 0. In other words, we shall consider wave- 
packets extended in z-direction, i.e. wave-ribbons. The 
eigenmodes are classified into transverse magnetic (TM) 
modes and transverse electric (TE) modes. For the TM 
modes, the magnetic field is in the xy-plane, while for the 
TE modes it is in the z-direction. The photonic bands 
are shown in Fig. E^b). The BC is along the z-direction 
in 2D, and its z-component is shown in Fig. It is 
seen that f2fc is strongly enhanced near the corners of the 



Brillouin zone, which can be interpreted as a 2D cut of 
the monopole structure in the extended space including 
the parameters, e.g. £ in the present case. [£ represents 
the degree of inversion-symmetry breaking.] 



TE 2nd band 




FIG. 4: Example of the BC of the 2D photonic crystal. It is 
prominent at the corners of the Brillouin zone. 

To induce the Hall effect of light in photonic crys- 
tals, we have to introduce a gradient of the refractive 
index, as seen in Eq. (J2J. Therefore, we shall intro- 
duce a slow variation of the envelope of the refractive 
index as n(r) — > n(r)/"/(x) where 1/7(2;) changes from 
1 to 1.2 within the range of w — 10a. We have solved 
the EOM for the wave-packets in the photonic crystal, 
r c = l{xc)^k a E nkc + k c x f2„ fec , k c = -[Vj(x c )]E nkc , 
where E nk and Q n k are the energy and the BC of the 
n-th band in the unperturbed crystal. The EOM for |z) 
brings about the phase shift for each mode. We set the 
initial k c near the corners of the Brillouin zone. This 
is because the BC is large at the corners and the effect 
by the anomalous velocity is expected to be prominent. 
The obtained trajectories are shown in Fig. EJ It is found 
that the shift of r c reaches to dozens of times the lattice 
constant especially for the TE second band. 



(a) Trajectory of center of gravity 




FIG. 5: (a) Real-space trajectory of the center of the wave- 
packet, (b) Momentum-space trajectory of the center of the 
wave-packet. The arrows in the same colors in (a) and (b) 
refer to the same photonic modes. 

We note that the BC is generally present in pho- 
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tonic crystals if the inversion-symmetry is broken. Our 
EOM offers an easy method to calculate the trajectory 
of wave-packets for generic photonic crystals By 
contriving crystal structures, this effect can be enhanced 
considerably as follows. The BC around the zone cor- 
ners is determined mostly by the splitting, 2|A|, be- 
tween the neighboring bands. Suppose, for example, 
at k = ko another band comes very close in energy to 
the one considered. The BC is approximately given by 
tl z ~ (A2+i,2|fc-fc |2)3/2 > where v is a nominal velocity 
of light. Thus when the light traverses this fco point, the 
shift is given by Sy c ~ sga\y Xc 7(x c )]v/ A. Therefore \Sy c \ 
is larger for smaller |A| as shown in Fig. El the second 
gap is smaller and hence the wave-packet in the second 
band shows larger \Sy c /a\. 

Let us discuss an effect of disorder in real photonic 
crystals. It has been discussed that the photonic band 
with lower index is more robust against disorder |2d| . 
Hence, in order to see the Hall effect of light, it is bet- 
ter to use the first band than the second band. While 
in Fig. [5] the shift of the first band is smaller than the 
second, it can be enhanced by reducing the gap between 
the two bands. The first band is not affected apprecia- 
bly when fluctuations of system parameters are less than 
a few percent. For visible or near infrared light, this 
requires fabrication precision within ~ 50nm, which is 
comparable to the current state of the art. Experimen- 
tally, the band splitting A of the order of 10 _5 nm _1 can 
be controlled in the photonic crystal for the visible light 
[2lj |. This corresponds to Sy c ~ 1/A ~ 100a. Thus, 
we can expect the shift as large as a hundred times the 
lattice constant a. 

To conclude, we have derived the semiclassical equa- 
tion of motion for the wave-packet of light including the 
Berry curvature and resultant anomalous velocity. This 
gives a natural generalization of geometrical optics in- 
cluding the wave nature of light, and leads to the Hall 
effect of light. This also offers the natural framework 
to describe the ray of light in the photonic crystal with- 
out the inversion symmetry where the Hall effect can be 
magnified tremendously. 
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Note added : A related argument on a spin transport of 
photons came to our attention after the submission of this 
work. K. Yu. Bliokh and Yu. P. Bliokh, |phy"sics/040211Ci| 
(unpublished) . 
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